Modelling Aotearoa New Zealand's COVID-19 protection framework and the transition away from the elimination strategy.

For the first 18 months of the COVID-19 pandemic, New Zealand used an elimination strategy to suppress community transmission of SARS-CoV-2 to zero or very low levels. In late 2021, high vaccine coverage enabled the country to transition away from the elimination strategy to a mitigation strategy. However, given negligible levels of immunity from prior infection, this required careful planning and an effective public health response to avoid uncontrolled outbreaks and unmanageable health impacts. Here, we develop an age-structured model for the Delta variant of SARS-CoV-2 including the effects of vaccination, case isolation, contact tracing, border controls and population-wide control measures. We use this model to investigate how epidemic trajectories may respond to different control strategies, and to explore trade-offs between restrictions in the community and restrictions at the border. We find that a low case tolerance strategy, with a quick change to stricter public health measures in response to increasing cases, reduced the health burden by a factor of three relative to a high tolerance strategy, but almost tripled the time spent in national lockdowns. Increasing the number of border arrivals was found to have a negligible effect on health burden once high vaccination rates were achieved and community transmission was widespread.

Using mechanistic model-based inference to understand and project epidemic dynamics with time-varying contact and vaccination rates.

Epidemiological models range in complexity from relatively simple statistical models that make minimal assumptions about the variables driving epidemic dynamics to more mechanistic models that include effects such as vaccine-derived and infection-derived immunity, population structure and heterogeneity. The former are often fitted to data in real-time and used for short-term forecasting, while the latter are more suitable for comparing longer-term scenarios under differing assumptions about control measures or other factors. Here, we present a mechanistic model of intermediate complexity that can be fitted to data in real-time but is also suitable for investigating longer-term dynamics. Our approach provides a bridge between primarily empirical approaches to forecasting and assumption-driven scenario models. The model was developed as a policy advice tool for New Zealand's 2021 outbreak of the Delta variant of SARS-CoV-2 and includes the effects of age structure, non-pharmaceutical interventions, and the ongoing vaccine rollout occurring during the time period studied. We use an approximate Bayesian computation approach to infer the time-varying transmission coefficient from real-time data on reported cases. We then compare projections of the model with future, out-of-sample data. We find that this approach produces a good fit with in-sample data and reasonable forward projections given the inherent limitations of predicting epidemic dynamics during periods of rapidly changing policy and behaviour. Results from the model helped inform the New Zealand Government's policy response throughout the outbreak.

Nanowire melting modes during the solid-liquid phase transition: theory and molecular dynamics simulations.

Molecular dynamics simulations have shown that after initial surface melting, nanowires can melt via two mechanisms: an interface front moves towards the wire centre; the growth of instabilities at the interface can cause the solid to pinch-off and breakup. By perturbing a capillary fluctuation model describing the interface kinetics, we show when each mechanism is preferred and compare the results to molecular dynamics simulation. A Plateau-Rayleigh-type of instability is found and suggests longer nanowires will melt via an instability mechanism, whereas in shorter nanowires the melting front will move closer to the centre before the solid pinch-off can initiate. Simulations support this theory; preferred modes that destabilise the interface are proportional to the wire length, with longer nanowires preferring to pinch-off and melt; shorter wires have a more stable interface close to their melting temperature, and prefer to melt via an interface front that moves towards the wire centre.

Real-time estimation of the effective reproduction number of SARS-CoV-2 in Aotearoa New Zealand.

During an epidemic, real-time estimation of the effective reproduction number supports decision makers to introduce timely and effective public health measures. We estimate the time-varying effective reproduction number, Rt , during Aotearoa New Zealand's August 2021 outbreak of the Delta variant of SARS-CoV-2, by fitting the publicly available EpiNow2 model to New Zealand case data. While we do not explicitly model non-pharmaceutical interventions or vaccination coverage, these two factors were the leading drivers of variation in transmission in this period and we describe how changes in these factors coincided with changes in Rt . Alert Level 4, New Zealand's most stringent restriction setting which includes stay-at-home measures, was initially effective at reducing the median Rt to 0.6 (90% CrI 0.4, 0.8) on 29 August 2021. As New Zealand eased certain restrictions and switched from an elimination strategy to a suppression strategy, Rt subsequently increased to a median 1.3 (1.2, 1.4). Increasing vaccination coverage along with regional restrictions were eventually sufficient to reduce Rt below 1. The outbreak peaked at an estimated 198 (172, 229) new infected cases on 10 November, after which cases declined until January 2022. We continue to update Rt estimates in real time as new case data become available to inform New Zealand's ongoing pandemic response.

Sensitivity of Reverse Transcription Polymerase Chain Reaction Tests for Severe Acute Respiratory Syndrome Coronavirus 2 Through Time.

BACKGROUND: Reverse transcription polymerase chain reaction (RT-PCR) tests are the gold standard for detecting recent infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Reverse transcription PCR sensitivity varies over the course of an individual's infection, related to changes in viral load. Differences in testing methods, and individual-level variables such as age, may also affect sensitivity. METHODS: Using data from New Zealand, we estimate the time-varying sensitivity of SARS-CoV-2 RT-PCR under varying temporal, biological, and demographic factors. RESULTS: Sensitivity peaks 4-5 days postinfection at 92.7% (91.4%-94.0%) and remains over 88% between 5 and 14 days postinfection. After the peak, sensitivity declined more rapidly in vaccinated cases compared with unvaccinated, females compared with males, those aged under 40 compared with over 40s, and Pacific peoples compared with other ethnicities. CONCLUSIONS: Reverse transcription PCR remains a sensitive technique and has been an effective tool in New Zealand's border and postborder measures to control coronavirus disease 2019. Our results inform model parameters and decisions concerning routine testing frequency.

An assessment of the potential impact of the Omicron variant of SARS-CoV-2 in Aotearoa New Zealand.

New Zealand delayed the introduction of the Omicron variant of SARS-CoV-2 into the community by the continued use of strict border controls through to January 2022. This allowed time for vaccination rates to increase and the roll out of third doses of the vaccine (boosters) to begin. It also meant more data on the characteristics of Omicron became available prior to the first cases of community transmission. Here we present a mathematical model of an Omicron epidemic, incorporating the effects of the booster roll out and waning of vaccine-induced immunity, and based on estimates of vaccine effectiveness and disease severity from international data. The model considers differing levels of immunity against infection, severe illness and death, and ignores waning of infection-induced immunity. This model was used to provide an assessment of the potential impact of an Omicron wave in the New Zealand population, which helped inform government preparedness and response. At the time the modelling was carried out, the date of introduction of Omicron into the New Zealand community was unknown. We therefore simulated outbreaks with different start dates, as well as investigating different levels of booster uptake. We found that an outbreak starting on 1 February or 1 March led to a lower health burden than an outbreak starting on 1 January because of increased booster coverage, particularly in older age groups. We also found that outbreaks starting later in the year led to worse health outcomes than an outbreak starting on 1 March. This is because waning immunity in older groups started to outweigh the increased protection from higher booster coverage in younger groups. For an outbreak starting on 1 February and with high booster uptake, the number of occupied hospital beds in the model peaked between 800 and 3,300 depending on assumed transmission rates. We conclude that combining an accelerated booster programme with public health measures to flatten the curve are key to avoid overwhelming the healthcare system.

Potential reduction in transmission of COVID-19 by digital contact tracing systems: a modelling study.

BACKGROUND: Digital tools are being developed to support contact tracing as part of the global effort to control the spread of COVID-19. These include smartphone apps, Bluetooth-based proximity detection, location tracking and automatic exposure notification features. Evidence on the effectiveness of alternative approaches to digital contact tracing is so far limited. METHODS: We use an age-structured branching process model of the transmission of COVID-19 in different settings to estimate the potential of manual contact tracing and digital tracing systems to help control the epidemic. We investigate the effect of the uptake rate and proportion of contacts recorded by the digital system on key model outputs: the effective reproduction number, the mean outbreak size after 30 days and the probability of elimination. RESULTS: Effective manual contact tracing can reduce the effective reproduction number from 2.4 to around 1.5. The addition of a digital tracing system with a high uptake rate over 75% could further reduce the effective reproduction number to around 1.1. Fully automated digital tracing without manual contact tracing is predicted to be much less effective. CONCLUSIONS: For digital tracing systems to make a significant contribution to the control of COVID-19, they need be designed in close conjunction with public health agencies to support and complement manual contact tracing by trained professionals.

A COVID-19 vaccination model for Aotearoa New Zealand.

We develop a mathematical model to estimate the effect of New Zealand's vaccine rollout on the potential spread and health impacts of COVID-19. The main purpose of this study is to provide a basis for policy advice on border restrictions and control measures in response to outbreaks that may occur during the vaccination roll-out. The model can be used to estimate the theoretical population immunity threshold, which represents a point in the vaccine rollout at which border restrictions and other controls could be removed and only small, occasional outbreaks would take place. We find that, with a basic reproduction number of 6, approximately representing the Delta variant of SARS-CoV-2, and under baseline vaccine effectiveness assumptions, reaching the population immunity threshold would require close to 100% of the total population to be vaccinated. Since this coverage is not likely to be achievable in practice, relaxing controls completely would risk serious health impacts. However, the higher vaccine coverage is, the more collective protection the population has against adverse health outcomes from COVID-19, and the easier it will become to control outbreaks. There remains considerable uncertainty in model outputs, in part because of the potential for the evolution of new variants. If new variants arise that are more transmissible or vaccine resistant, an increase in vaccine coverage will be needed to provide the same level of protection.

Effect of vaccination, border testing, and quarantine requirements on the risk of COVID-19 in New Zealand: A modelling study.

We couple a simple model of quarantine and testing strategies for international travellers with a model for transmission of SARS-CoV-2 in a partly vaccinated population. We use this model to estimate the risk of an infectious traveller causing a community outbreak under various border control strategies and different levels of vaccine coverage in the population. Results are calculated from N = 100,000 independent realisations of the stochastic model. We find that strategies that rely on home isolation are significantly higher risk than the current mandatory 14-day stay in government-managed isolation. Nevertheless, combinations of testing and home isolation can still reduce the risk of a community outbreak to around one outbreak per 100 infected travellers. We also find that, under some circumstances, using daily lateral flow tests or a combination of lateral flow tests and polymerase chain reaction (PCR) tests can reduce risk to a comparable or lower level than using PCR tests alone. Combined with controls on the number of travellers from countries with high prevalence of COVID-19, our results allow different options for managing the risk of COVID-19 at the border to be compared. This can be used to inform strategies for relaxing border controls in a phased way, while limiting the risk of community outbreaks as vaccine coverage increases.

Early intervention is the key to success in COVID-19 control.

New Zealand responded to the COVID-19 pandemic with a combination of border restrictions and an Alert Level (AL) system that included strict stay-at-home orders. These interventions were successful in containing an outbreak and ultimately eliminating community transmission of COVID-19 in June 2020. The timing of interventions is crucial to their success. Delaying interventions may reduce their effectiveness and mean that they need to be maintained for a longer period. We use a stochastic branching process model of COVID-19 transmission and control to simulate the epidemic trajectory in New Zealand's March-April 2020 outbreak and the effect of its interventions. We calculate key measures, including the number of reported cases and deaths, and the probability of elimination within a specified time frame. By comparing these measures under alternative timings of interventions, we show that changing the timing of AL4 (the strictest level of restrictions) has a far greater impact than the timing of border measures. Delaying AL4 restrictions results in considerably worse outcomes. Implementing border measures alone, without AL4 restrictions, is insufficient to control the outbreak. We conclude that the early introduction of stay-at-home orders was crucial in reducing the number of cases and deaths, enabling elimination.

Vaccination and testing of the border workforce for COVID-19 and risk of community outbreaks: a modelling study.

Throughout 2020 and the first part of 2021, Australia and New Zealand have followed a COVID-19 elimination strategy. Both countries require overseas arrivals to quarantine in government-managed facilities at the border. In both countries, community outbreaks of COVID-19 have been started via infection of a border worker. This workforce is rightly being prioritized for vaccination. However, although vaccines are highly effective in preventing disease, their effectiveness in preventing infection with and transmission of SARS-CoV-2 is less certain. There is a danger that vaccination could prevent symptoms of COVID-19 but not prevent transmission. Here, we use a stochastic model of SARS-CoV-2 transmission and testing to investigate the effect that vaccination of border workers has on the risk of an outbreak in an unvaccinated community. We simulate the model starting with a single infected border worker and measure the number of people who are infected before the first case is detected by testing. We show that if a vaccine reduces transmission by 50%, vaccination of border workers increases the risk of a major outbreak from around 7% per seed case to around 9% per seed case. The lower the vaccine effectiveness against transmission, the higher the risk. The increase in risk as a result of vaccination can be mitigated by increasing the frequency of routine testing for high-exposure vaccinated groups.

Structure of the Region-Technology Network as a Driver for Technological Innovation.

Agglomeration and spillovers are key phenomena of technological innovation, driving regional economic growth. Here, we investigate these phenomena through technological outputs of over 4,000 regions spanning 42 countries, by analyzing more than 30 years of patent data (approximately 2.7 million patents) from the European Patent Office. We construct a bipartite network-based on revealed comparative advantage-linking geographic regions with areas of technology and compare its properties to those of artificial networks using a series of randomization strategies, to uncover the patterns of regional diversity and technological ubiquity. Our results show that the technological outputs of regions create nested patterns similar to those of ecological networks. These patterns suggest that regions need to dominate various technologies first (those allegedly less sophisticated), creating a diverse knowledge base, before subsequently developing less ubiquitous (and perhaps more sophisticated) technologies as a consequence of complementary knowledge that facilitates innovation. Finally, we create a map-the Patent Space Network-showing the interactions between technologies according to their regional presence. This network reveals how technology across industries co-appear to form several explicit clusters, which may aid future works on predicting technological innovation due to agglomeration and spillovers.

Maori and Pacific people in New Zealand have a higher risk of hospitalisation for COVID-19.

AIMS: We aim to quantify differences in clinical outcomes from COVID-19 infection in Aotearoa New Zealand by ethnicity and with a focus on risk of hospitalisation. METHODS: We used data on age, ethnicity, deprivation index, pre-existing health conditions and clinical outcomes on 1,829 COVID-19 cases reported in New Zealand. We used a logistic regression model to calculate odds ratios for the risk of hospitalisation by ethnicity. We also considered length of hospital stay and risk of fatality. RESULTS: After controlling for age and pre-existing conditions, we found that Maori have 2.50 times greater odds of hospitalisation (95% CI 1.39-4.51) than non-Maori non-Pacific people. Pacific people have three times greater odds (95% CI 1.75-5.33). CONCLUSIONS: Structural inequities and systemic racism in the healthcare system mean that Maori and Pacific communities face a much greater health burden from COVID-19. Older people and those with pre-existing health conditions are also at greater risk. This should inform future policy decisions including prioritising groups for vaccination.

Managing the risk of a COVID-19 outbreak from border arrivals.

In an attempt to maintain the elimination of COVID-19 in New Zealand, all international arrivals are required to spend 14 days in government-managed quarantine and to return a negative test result before being released. We model the testing, isolation and transmission of COVID-19 within quarantine facilities to estimate the risk of community outbreaks being seeded at the border. We use a simple branching process model for COVID-19 transmission that includes a time-dependent probability of a false-negative test result. We show that the combination of 14-day quarantine with two tests is highly effective in preventing an infectious case entering the community, provided there is no transmission within quarantine facilities. Shorter quarantine periods, or reliance on testing only with no quarantine, substantially increases the risk of an infectious case being released. We calculate the fraction of cases detected in the second week of their two-week stay and show that this may be a useful indicator of the likelihood of transmission occurring within quarantine facilities. Frontline staff working at the border risk exposure to infected individuals and this has the potential to lead to a community outbreak. We use the model to test surveillance strategies and evaluate the likely size of the outbreak at the time it is first detected. We conclude with some recommendations for managing the risk of potential future outbreaks originating from the border.

Estimated inequities in COVID-19 infection fatality rates by ethnicity for Aotearoa New Zealand.

AIMS: There is limited evidence as to how clinical outcomes of COVID-19 including fatality rates may vary by ethnicity. We aim to estimate inequities in infection fatality rates (IFR) in New Zealand by ethnicity. METHODS: We combine existing demographic and health data for ethnic groups in New Zealand with international data on COVID-19 IFR for different age groups. We adjust age-specific IFRs for differences in unmet healthcare need, and comorbidities by ethnicity. We also adjust for life expectancy reflecting evidence that COVID-19 amplifies the existing mortality risk of different groups. RESULTS: The IFR for Maori is estimated to be 50% higher than that of non-Maori, and could be even higher depending on the relative contributions of age and underlying health conditions to mortality risk. CONCLUSIONS: There are likely to be significant inequities in the health burden from COVID-19 in New Zealand by ethnicity. These will be exacerbated by racism within the healthcare system and other inequities not reflected in official data. Highest risk communities include those with elderly populations, and Maori and Pacific communities. These factors should be included in future disease incidence and impact modelling.

Stability of amphiphilic Janus dimers in shear flow: a molecular dynamics study.

Amphiphilic Janus particles in a flow are thought to experience a torque due to the asymmetry in slip at their surfaces. This effect has the potential to destabilise self-assembled Janus structures in flows due to the forces and torques applied to individual Janus nanoparticles. In this work, we investigate the stability of amphiphilic Janus dimers and homogeneous hydrophobic dimers in shear flow using molecular dynamics, and study possible break-up mechanisms. In particular, we consider the influence of the activation enthalpy and entropy on the thermal break-up rate of these dimers. Janus dimers are less stable than hydrophobic dimers, and increasing the applied shear rate has a greater effect on break-up for Janus dimers. Two mechanisms leading to increased break-up in shear flow are studied, namely the rotational speed of the dimers and the orientation of individual spheres in the dimers, and we propose a descriptive equation for calculation of the break-up rate. Overall, the results indicate that the stability of dimers in shear flow depends on the slip length at the spheres' surfaces, and that the slip length difference on Janus dimers could contribute to destabilisation.

Molecular dynamics simulations of Janus nanoparticles in a fluid flow.

We study the forces and torques on individual Janus nanoparticles in a fluid flow using molecular dynamics simulations. In particular, we consider amphiphilic Janus nanospheres that have different slip boundary conditions on each hemisphere, and calculate the forces and torques experienced by them as a function of their orientation with respect to the flow direction. Furthermore, we examine nanoparticles that are deformed slightly from a spherical shape, and have no-slip boundary conditions. We compare the simulation results to the predictions of previously introduced theoretical approaches, which compute the forces and torques on particles with variable slip lengths or aspherical deformations that are much smaller than the particle radius. We found that there is good qualitative agreement between the forces and torques computed from our simulations and the theoretical predictions, and the forces quantitatively agree when the assumptions made in the theoretical descriptions are satisfied.

Surface melting and breakup of metal nanowires: Theory and molecular dynamics simulation.

We consider the surface melting of metal nanowires by solving a phenomenological two-parabola Landau model and by conducting molecular dynamics simulations of nickel and aluminum nanowires. The model suggests that surface melting will precede bulk melting when the melt completely wets the surface and the wire is sufficiently thick, as is the case for planar surfaces and sufficiently large nanoparticles. Surface melting does not occur if the melt partially wets or does not wet the surface. We test this model, which assumes that the surface energies of the wire are isotropic, using molecular dynamics simulations. For nickel, we observe the onset of anisotropic surface melting associated with each of the two surface facets present, but this gives way to uniform surface melting and the solid melts radially until the solid core eventually breaks up. For aluminum, while we observe complete surface melting of one facet, the lowest energy surface remains partially dry even up to the point where the melt completely penetrates the solid core.

Control of nanoparticle formation using the constrained dewetting of polymer brushes.

We have used coarse-grained molecular dynamics simulations to investigate the use of pinned micelles formed by the constrained dewetting of polymer brushes to act as a template for nanoparticle formation. The evaporation of a thin film containing a dissolved solute from a polymer brush was modeled to study the effect of solubility, concentration, grafting density, and evaporation rate on the nucleation and growth of nanoparticles. Control over particle nucleation could be imposed when the solution was dilute enough such that particle nucleation occurred following the onset of constrained dewetting. We predict that nanoparticles with sizes on the order of 1 nm to 10 nm could be produced from a range of organic molecules under experimentally accessible conditions. This method could allow the functionality of organic materials to potentially be imparted onto surfaces without the need for synthetic modification of the functional molecule, and with control over particle size and aggregation, for application in the preparation of surfaces with useful optical, pharmaceutical, or electronic properties.